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We consider the dynamics of the charge carriers in a tunneling- enhanced percolation network, 
named as a Random Resistor cum Tunneling-bond Network (RRTN), where we allow tunneling in 
the gap between two randomly thrown nearest neighbour metallic bonds only. Our earlier stud- 
ies involved the dc and ac nonlinear response, the percolative aspects, dielectric breakdown, low- 
temperature variable range hopping (VRH) conduction, etc. in the RRTN. Here we study the 
non-equilibrium dynamics of the carriers. With two far-from-equilibrium, initial inverse power-law 
relaxations extending over several decades, the dynamics has a lot of similarities with a wide va- 
riety of naturally occuring avalanche-like, run-away phenomena in driven, disordered systems with 
statistically correlated randomness. In the power-law regime, the RRTN violates the Boltzmann's 
(or, Debye) relaxation time approximation strongly. Beyond this regime, the response relaxes ex- 
ponentially fast (acquires one time-scale) to a steady-state, and thus the relaxation approximation 
becomes exact. 

PACS numbers: 71.55.Jv, 89.75.-k, 89. 75. He 



I, INTRODUCTION 



Slow relaxation phenomena in random composite materials made of components having widely different generalized 
susceptibilities {e.g., permeability, dielectric constant, electrical/ thermal conductivity, viscosity, elastic modulii etc.), 
continue to remain intriguing and hence a topic of intense research. The susceptibility, being a measure of the 
response of a system to an appropriate external perturbation or a driving force, is typically assumed to be linear 
under a vanishingly small external force. But, many natural systems/ phenomena do not show any measurable 
response until an appropriate driving force exceeds a measurable threshold. For example, a rigid body on a rough 
surface does not move until the driving force exceeds a finite frictional force. More intriguingly, it is almost a law 
of nature that in such driven, macroscopic systems, the response characteristic is nonlinear Q with a concommitant 
critical behaviour at the threshold. Besides the case of frictional motion, nonlinear response above a threshold is also 
found in pinned charge-density- wave systems, disordered granular superconductors Q, etc. 

In studying the dynamical behaviour of a material, one usually measures its appropriate response property, say, 
4>{t), as a function of time t, from a non-equilibrium to an equilibrium (for a closed system) or to a steady (for an open 
system) state. In general, this relaxation is classified into two groups: (i) a purely Debye type with an exponential 
relaxation function, 4>{t) = exp(— t/r), r being a characteristic time-scale, called the relaxation time; or (ii) a non- 
Debye type where 4>(t) is either a linear superposition of exponential functions, or a sub-exponential function (as in 
some glassy systems), with multiple relaxation times. In more intriguing cases, this non-Debye relaxation may even 
be a power-law or a logarithmic function, without any time-scale, or one with r — > oo. The appearence of such a 
scale-free, slow dynamics with one or more power-laws from the early stages of evolution, is what concerns us here 
because of the fundamental issue of the failure of Boltzmann's relaxation time approximation Q in driven systems 
far-from-equilibrium. 



II. EXPERIMENTS AND SOME RELATED MODELS 

In a lucidly written review Scher, Shlesinger and Bendler 0| focus on experimental observations (1970's onwards) 
and the origin of two power-law kinetics. For a few examples, we cite some transient photocurrent measurements Q 
on a-Si:H, a-As2Se3 etc. and a couple of hole transport t 6j data on PVK and Si-MOS devices, where two consecutive 
power-law decays, of the forms t~ a and t~® (a, (3 > 0), covering one or more decades in time each (with a crossover in- 
between) were observed. Based on the continuous time random walk with a long-tailed power-law probability density 
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(a) o-bonds in a p=0.30, L=20 RRN sample 
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(b) RRTN based on the RRN beside: o-bond 
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Figure 1: (a) A typical configuration of random metallic o-bonds in a L = 20, p = 0.3 RRN sample on a square lattice. 
The insulating bonds are kept blank. This RRN configuration, being below its threshold of p c = 0.5, is not percolating (an 
insulating RRN); (b) In the ensuing maximal (with fully statistically correlated tunneling t-bonds) RRTN configuration, the 
t-bonds (dashed lines) are placed between any two nearest neighbour o-bonds. They get activated when the potential differences 
across them exceed a microscopic voltage threshold, v g . With all the allowed t-bonds active, this maximal RRTN percolates 
(is a 'metallic' RRTN). 



function for the random waiting times (release time of trapped carriers by tunneling) , Scher et al. [J] formulated a 
theory regarding the above results. The latter long-tailed power-law function violates the Central Limit Theorem, 
since all of its moments including the first (mean waiting-time) diverge. The unifying feature of the above random 
walk is the scale-invariance of the shape of the relaxation current I(t), if one normalizes the time by a transit time 
t r , which is a sample dependent parameter. This stochastic theory by Scher et al. |4] explains the results of many 
early experiments, following the relation a + (3 = 2. 

But, there is a huge variety of relatively recent, more intriguing experiments on soft-condensed or complex systems 
where the couple of exponents a and (3 do not seem to follow any simple algebraic relation. Weron and Jurlewicz note 
that an experiment on dielectric relaxation in a system of dipoles 7] involves the crossover between different forms 
of self-similarity. They argue that 0] a couple of power-law decays appear due to the coupling of micro-clusters of 
dipoles with a distribution of r's (i.e., multiple time-scales). An experiment on fluorescence intermittency of single 
ZnS overcoated CdSe quantum dots the distribution of on and off times (blinking kinetics) is reported to follow a 
single power-law [2^|. In a biological system of photo-dissociated heme-proteins, the rebinding of the ligands of iron 
(i.e., the CO and the O2 molecules) is observed to follow an inverse power-law dynamics [10]. In a theoretical study 
of the same work, Tsallis et al. [Tl| claim two inverse power-law regimes and demonstrated its possible connection 
with nonextensive thermo- statistics (entropy). Naudts and Czachor ^3] analyzed the data of many experiments 
including those of 0, 0] , and maintain that these two power-law decays result from some choice of parameters of 
their probability density function. All of the above theories lead to the result that a < [3, if two power-law relaxations 
exist. In a different biological system, the dynamics of Ca channels in living cells 0], the distribution of the survival- 
times of the channels has been studied. A stochastic dynamical model, with one dimensional geometry, was proposed 
specifically to explain only the later power-law dynamics (2r| . 

Two power-law growths are also claimed to have been found, using atomic force microscopy (AFM), in the early 
dynamics of sputtered Ag particles on Si(0 1) substrate These authors relate the two different growth-dynamics 
to two competing structural rearrangements at two different length-scales. In the case of the growth of a large single 
DLA (diffusion limited aggregation) cluster, using upto 10 s particles in a computer experiment, two power-laws seem 
to dictate the growth process as a function l6j of a time-like entity. A similar theme has been reflected in various 
methods (e.g., the sol-gel method) of synthesis of nanomaterials This involves a slow relaxation of restructuring 
many local clusters by crossing the local barriers due to a caging effect (typical of liquid-like, amorphous or dense 
granular materials) followed by (or, may be, added to) another slow relaxation of global restructuring of the local 
clusters themselves. Even in some relatively recent studies of cellular automata models of earthquakes [la j. one 
power-law dynamics seems to result under certain choice of parameters. 
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III. A TUNNELING PERCOLATION MODEL 

To study the basic physics behind the transport properties within a minimal model, Sen and co-workers [TflL l20T | 
constructed a lattice based bond percolation model. In this model, in addition to the randomly placed ohmic bonds 
(linear response) in the Random Resistor Network (RRN), tunneling bonds (t-bond) were introduced between two 
nearest neighbour (nn) ohmic bonds (o-bonds) separated by one insulator only. This tunneling is actually semi- 
quantum or semi-classical in the sense that no quantum mechanical phase-information of the charge carriers appears 
in the model. In the steady-state, a t-bond is insulating if the microscopic voltage difference across it (v) is less than 
a fixed voltage threshold (v g , identical for all the t-bonds) and conducts linearly a current if v > v g . Indeed, the 
phenomenological parameters p and v g ensure that both the disorder and the coulomb interaction (nonlinear response 
[? ] as an outcome) is in-built in this model. 

Historically 19], we call this tunneling-enhanced percolation model as the Random Resistor cum Tunneling Network 
(RRTN). The appearence of the t-bonds in this perfectly correlated, i.e., deterministic fashion is the origin of a very 
low percolation threshold, in the RRTN model. For example, in the fig^a), the RRN with p=0.3, does not percolate. 
In this respect, one may note that the percolation threshold for the RRN (square lattice) is p c = 0.50. But, the 
maximal RRTN (where all the possible t-bonds, both along and across the field, are active, i.e., conducts a maximum 
possible current in the steady state) generated from that particular RRN, percolates as shown in the figQIb). Standard 
calculations with finite-size scaling analysis, establishes a new percolation threshold of p c t = 0.18 j^, for the maximal 
RRTN. 

Further, because of the finite threshold of the t-bonds, the bulk dc conductance G(V) (defined either as y or as 
jjy) is a strongly nonlinear curve, as stated earlier. At a sufficiently low voltage (below v g ), no t-bond is active, 
and hence the effective bulk conductance is only due to the percolating ohmic backbone. So the G(V) vs V graph 
is parallel to the V-axis (falls on the T^-axis if the ohmic backbone is non-percolating). This region is termed as 
the lower linear/ohmic regime. Thereafter, under further increment of the external voltage V, yet more microscopic 
t-bonds cross the threshold (v g ), i.e., become active, and give rise to extra parallel paths in the network. At an very 
high voltage [2lJ (asymptotically infinite), all the possible t-bonds become active (i.e., creates a maximal RRTN) and 
hence the conductance saturates to its largest possible value (for that particular configuration) . Obviously, this region 
belongs to the upper linear/ohmic regime. Other than the success of this discrete model in understanding various 
experiments on an ultra-low percolation threshold, and various aspects of nonlinear dc and ac responses, it has also 
been quite useful for understanding some very unusual aspects observed in low-temperature variable range hopping 
(VRH) conduction [22], and some interesting aspects of breakdown phenomena (23 . Thus, even though, time enters 
in an implicit fashion in some of the above studies, an explicit characterization of the relaxation dynamics in the 
RRTN model in the perspective of various experiments and models (as described above) was considered necessary. 

IV. NON-EQUILIBRIUM DYNAMICS IN THE RRTN MODEL 

In this paper, we study the current relaxation dynamics towards a steady-state in the RRTN model. Away from 
equilibrium, a t-bond with a microscopic v < v g behaves like a dielectric material between two metals (o-bonds); and 
the resulting charging effect gives rise to a displacement current (C^&-, where C is the capacitance). Thus, a t-bond 
gives rise to a displacement current if v(t) < v g , and an 'ohmic' current if v(t) > v g j24[. For our calculations, we use 
the values of the microscopic conductance g = 1.0 (o-bonds), gt — 10 -2 (t-bonds), v g = 0.5, and C — 10~ 5 for the 
t-bonds (in some arbitrary units). 

In our numerical study, we apply an uniform electric field across RRTN's of different system sizes (L) and ohmic 
bond concentrations (p) . We study the evolution of the current in a RRTN starting from the switching on state until 
it approaches its asymptotic steady-state. To do this, we follow basically the current conservation (Kirchhoff's laws) 
locally at each node of the lattice. The aim is to study the achievement of a global current conservation as an outcome 
of the local current conservation (hence, the dynamics). A discrete, scaled time unit has been chosen as one complete 
scan through each site of the lattice. This local conservation or the equation of continuity reads as, 

X>„(*) = 0, Vi. (1) 

Here the sum has been taken over currents i nn (t) through various types of nearest neighbour (nn) microscopic bonds 
around any node/site of the lattice. For the case of a square lattice, one considers the four m's around a node 
inside the bulk (three and two nn's respectively at any boundary or a corner). If eq. were true simultaneously for 
each site of the lattice, then the global conservation (the steady state) for the entire network would automatically be 
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Figure 2: Non-equilibrium dynamics in terms of a scaled time shows two initial power-law decays (with intermediate crossover) 
in the current relaxation for (a) one RRTN sample (p = 0.55, L — 60) with the exponents of 0.44 and 0.72; and (b) another 
sample (p = 0.35, L = 80) where the exponents are 0.78 and 0.57 respectively. The dynamics in the RRTN is always followed 
by an exponentially decaying tail (not shown in the figures here) towards its steady state. 

achieved. As we need to start with an initial (arbitrary) microscopic voltage distribution, the eq. would not hold 
for all the sites of the lattice. Some correction term would be required at each site and this requirement leads to the 
following time evolution algorithm which we call as the lattice Kirchhoff 's dynamics: 

vv,k,t+i)=vv,k,t)+m^$-, (2) 

/ ) 9nn 

where g nn are the various microscopic conductances of the nn bonds around the node (j, k) . Then we numerically 
solve a set of coupled difference equations on the lattice. The move towards a macroscopic steady-state implies that 
the difference of currents through the first and the last layers tends to zero as a function of time. In practice, the 
system is considered to have reached its steady state when this difference decreases to a pre-assigned smallness. 

V. RESULTS AND DISCUSSIONS 

As shown in the figs. El we observe a non-Debye type current relaxation with two consecutive initial power-laws 
(and a crossover in-between), each of them spanning at least a decade. This asymptotic steady state current (whether 
insulating or conducting) for any randomly chosen RRTN, is found to be very robust against any initiallly chosen 
voltage distribution on the lattice. To analyse the current evolution of an insulating or a conducting sample on the 
same footing, we subtract the corresponding steady current, I(t — > oo), from the evolving current I(t) at the time 
t. Our observation on the transient current response, indicates clearly the existence of a couple of initial power-laws, 
whose exponents differ significantly for systems with different configurations, p, L and external voltage V. So, we 
choose to work with one sample at a time and analyze its results within our numerical accuracy. 

For example, in the fig-Gfa), we show the dynamics for a sample with L = 60, p — 0.55, V = 2.0 and in the figE^b) 
another with L = 80, p = 0.35, V = 20.0. The first figure (i.e., fig-GIa)) represents the class of relaxation, where the 
second exp onent (0.72) is larger than the first (0.44) (the only class reported in our previous work at a particular 
p = p c Hif). There are quite a few theoretical works 0,0, EUl in this regard. Some other experiments (l^.ll4.IT5| find 
a second exponent smaller than the first. Since our earlier report we have been able to reproduce this other class 
of relaxation (e.g., fig-E^b) with the exponents 0.78, 0.57), as well, within the context of our RRTN. In special cases, 
we do find only one power-law relaxation, which may be considered to be the borderline between the above two, or as 
the merging of the two power-law exponents. Further, we do not find any particular relation between the exponents. 
Next, the existence or the lack of any asymptotic exponential kinetics is not explicitly stated in most of them. Indeed, 
in some of the theoretical studies (e.g. in |j,@,^|), the second power-law persists upto asymptotically infinite times. 
This trend can not describe the possiblity for these systems to reach an appropriate steady state. In contrast, beyond 
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the power-law relaxations (whether one or two) , our model acquires a relaxation time r and the system follows a fast 
exponential dynamics as a precursor of a final steady state which is effectively diffusive (or ohmic). Obviously, the 
power-law relaxations at times t < r imply a strong deviation from the Boltzmann's relaxation time approximation 
(i.e., strongly non-Debye type). 

As far as the ori gin of the two power-law dynamics are concerned, we have already outlined the main content of 
some of them P.l^.llllll2.ll3.ll^| in the section on experiments. In most of them they occur due to local structural 
rearrangements preceding the final global structural rearrangements. In our case, the basic structure in a particular 
sample is created once for all, and it is the fields across the bonds which keep changing in such a fashion that the local 
conservation (Kirchhoff 's laws at each node) dominates the first power-law regime and the global current conservation 
is dominated in the latter power-law regime (of course, there are the required structural rearrangements of the active 
t-bonds). So, it is interesting that these two very different mechanisms give rise to a qualitatively identical outcome. 
Further, since the power-law dynamics occurs even for p's away from p c or p c t, it clearly demonstrates that they are 
not organized by any type of criticality. Finally, as discussed above, while most of the other theoretical works, are 
destined to get only one class of two-power-law relaxation behaviour (namely, a < (3), the RRTN dynamics has the 
ability to capture both of the classes for different sets of parameters. 
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